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Abstract 

Near-critical binary mixtures containing ions and confined between two charged and selective 
surfaces are studied within a Landau-Ginzburg theory extended to include electrostatic interactions. 
Charge density profiles and the effective interactions between the confining surfaces are calculated 
in the case of chemical preference of ions for one of the solvent components. Close to the consolute 
point of the binary solvent, the preferential solubility of ions leads to the modification of the 
charge density profiles in respect to the ones obtained from the Debye-Hiickel theory. As a result, 
the electrostatic contribution to the effective potential between the charged surface can exhibit 
an attractive well. Our calculations are based on the approximation scheme valid if the bulk 
correlation length of a solvent is much larger than the Debye screening length; in this critical 
regime the effect of charge on the concentration profiles of the solvent is subdominant. Such 
conditions are met in the recent measurements of the effective forces acting between a substrate 
and a spherical colloidal particle immersed in the near-critical water- lutidine mixture [Nature 451, 
172 (2008)]. Our analytical results are in a quantitative agreement with the experimental ones. 

PACS numbers: 05.70.Np,05.70.Jk,82.45.Gj 
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I. INTRODUCTION 



Ions dissolved in a binary liquid mixture often display preferential solubility in the one 
component of the solvent. Also the interaction between the ± ions and the solvent can be 
different, which is known as an unequal partitioning of ions. In a bulk system, a selective 
solvation leads to the shift of the critical point of the demixing transition l|-|3| and to a 



number of other effects on the phase separation of a binary solvent 



^. In the presence of 



external charged surfaces the shift of the bulk critical point of a mixture can be enhanced 



by a dielectric inhomogeneity arising due to t 
the charged surface (dielectrophoretic forces) 



le attraction of high permittivity solvent to 
5|-[7|. Moreover, a selective solvation of ions 



may change the concentration profiles of the binary solvent near the wall and, reversely, 
adsorption phenomena can significantly influence the distribution of ions near a charged 
surface. These mutual influences have been recently studied theoretically for the case when 
the binary solvent is near its consolute point in the semi-infinite geometry 8|, and for systems 
confined between two parallel walls or substrates IJIlll. In the latter case, the consequences 
of the ions-solvent coupling for the effective forces acting on the confining surfaces were 
studied. 

One of the motivation for such investigations is provided by recent experimental works 



12h14|. where the effective potential between a charged substrate and a likely charged 



colloidal particle immersed in a water-lutidine critical mixture (T^ ^ 307. 15 J^) was directly 
measured. The surfaces of the colloidal particle and of the flat substrate with similar or 
opposite adsorption preferences were used in order to verify predictions of the theory for the 
thermodynamic Casimir force. These, so called, critical Casimir forces acting between the 
colloidal particle and a flat substrate arise as a result of the modiflcations of the relevant 
order parameter (OP) and restrictions of its fluctuation spectrum by the conflning surfaces. 
Close to the critical point of the solvent, attraction is predicted for like surfaces, whereas 
repulsion is predicted if one surface is hydrophilic and the other one is hydrophobic. 

The theory of effective interactions between two surfaces conflning a near-critical fluid 



is well developed for uncharged surfaces and for mixtures of neutral components 



15 



18|. 



However, in the experiments mentioned above the surfaces were charged, and moreover 



a small amount of ions was present in the solution. In Ref. 12, 131 the ions result from 

n 

dissociation of water, and in Ref. [1^ a hydrophilic salt was added. Far from Tc repulsion 



2 



has been present independently o: 
electrostatic potential dominates 



" the adsorption preferences of the surfaces, because the 
12Nl4|. The electrostatic repulsion decays exponentially 
with the decay rate equal to the Debye screening length 1/k. For T — )■ Tc, in addition to 
the repulsion for small separations L ~ an attraction (repulsion) has been observed 
for surfaces with the same (opposite) adsorption preferences for larger separations, L ~ ^, 



12 



13| . Such behavior is predicted 



where ^ is the bulk correlation length of the solvent 
by the sum of the electrostatic and the critical Casimir potentials for the corresponding 
boundary conditions. A full quantitative agreement between the experiment and the sum of 



the electrostatic and the critical Casimir potentials could not be obtained, however |12l. |13|. 
The sum of the electrostatic potential that fitted well the experimental results far from 
and of the critical Casimir potential that fitted well the data for separations L ^ 1/k, close to 
Tc, for intermediate distances disagreed strongly with the measured potential. The authors 
concluded that coupling between the critical concentration fluctuations and the distribution 
of ions may lead to modifications of the potential. For this reason, only distances significantly 
larger than the screening length were considered close to the critical point to verify the 
theory of the critical Casimir potential. In the presence of salt a more complex behavior 
was observed, in particular, an unexpected attraction between hydrophilic and hydrophobic 
surfaces for intermediate temperatures (^k < 1) UlnlQ. 

The shape of the effective interaction potentials between charged selective surfaces con- 
fining the critical binary solvent with ions resembles strongly the intermolecular interaction 
potentials, but on a much larger scale. Because of possible applications, the ability to design 
the interaction potential between, e.g., the two colloidal particles is of interest, therefore the 
mutual effect of ion distribution and concentration profiles deserves serious attention. 



Theoretical studies reported in Refs. |9l-lll 



19| are all based on the Ginzburg-Landau-like 



theory but with a different level of complexity as far as the parameter space is concerned. 
In general, a high-dimensional parameter space is required for a full description of a four- 
component mixture, with two of the species carrying a charge, in a presence of two charged 



and selective surfaces. In Refs. 



10 



11| the reduced description has been employed in order to 



investigate the particular mechanisms and the role of the specific interactions. Accordingly, 
in these studies, e.g., the van der Waals (vdW) type of interactions between ions and be- 



tween ions and the walls have been neglected altogether. In Ref. 



10| , a non-trivial interplay 



between critical and electrostatic phenomena (which goes beyond the simple superposition 
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of the critical Casimir and the electrostatic potentials) arises as a result of an unequal parti 



tioning of the salt ions in a non-uniform solvent. In Ref. the focus is on the electrostatic 
effects, therefore also interactions between the components of the solvent and the walls have 
been neglected. A preference of charged walls for one of the solvent components (with the 
largest permittivity) has been taken into account via the composition-dependent permittiv- 
ity. Within this approach, for an equal partitioning of the salt ions in each component of the 
solvent, an attraction between like-charge surfaces can occur as a result of dielectrophoretic 
forces and the ion-solvent coupling. 

Here we extend the theoretical approach developed in Ref. 8j for a semi-infinite system to 
the slit geometry and determine the influence of critical adsorption on the charge distribution 
close to the critical point of the solvent, i.e., for k,^ > 1. Such a ratio of relevant length 
scales in the system has been realized in the experiments described in Ref. 



12 



13|. Next 



we examine the effect of these modifications of the distribution of ions on the form of the 
effective potential between confining surfaces which are charged and selective. Within the 
approximation scheme that we use in our analysis, this effective potential can be written as a 
sum of three contributions: the critical Casimir potential, the pure electrostatic potential (as 
given by the linearized Debye-Hiickel (DH) theory), and the potential arising from the ion- 
solvent coupling. We use a Derjaguin approximation 20| in order to compare our theoretical 
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13|. 



predictions for the effective potential with the experimental data reported in Ref. 
Within the Derjaguin approximation the interaction potential between the sphere and the 
planar wall is expressed in terms of the interaction potential in the slit. For the critical 
Casimir part of the total effective potential we use the scaling function determined to a 
great degree of accuracy from the MC simulations in d = 3 21 1. 



The description of the system used in the present paper is more complete than the ones 
used in Refs. [lO, ll], 19 1 in the sense that it treats ions as the molecules which interact 



also non-electrostatically with each other and with the walls. Consequently, we consider a 
system confined by two charged walls which are selective to all components of the mixture. 
Due to the more complete description, the parameter space of our model is somewhat larger 

n n n n 

than in the other approaches [10|, llll, 119|. In the full version of the model |8[, the vdW 
interactions between all pairs of components of the mixture, and the dependence of the 
permittivity on the concentration was assumed. The number of parameters can be reduced 
for particular systems. For example, for hydrophilic ions we are left with 3 parameters 



4 
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characterizing non-Coulombic interactions [8|], while two such parameters are present in Refs. 



Rather general description developed from microscopic theory lends itself to still 



another mechanisms leading to the unintuitive effects in the slit geometry [9|. Moreover, 
the general framework of our theory is also suitable for antagonistic salt, which leads to 
interesting phenomena {22]. Finally, because the Ginzburg-Landau-type theory that we 
employ has been developed from the microscopic lattice gas model of the four- comp onent 



mixture, in our model the entropy of mixing is better approximated than in Refs. [10|, 
where the entropy of mixing has been taken separately for the binary solvent (without 
ions) and separately for the ions (as an entropy of an ideal gas). In the present work we 



assume, as in Ref. 



10| a uniform permittivity, because the dielectrophoretic effects can be 



mimicked by an appropriate contributions to the surface fields. Here we consider the case 
of the equal partitioning of the ions in the solvent (like in Ref. ll|). What distinguishes 
our study from the other similar approaches proposed recently, is that our analytical results 
are obtained beyond the linear approximation for the EL equations, and a quantitative, not 
only a qualitative agreement with experiments is obtained. As in Ref. [loj] , interesting effects 
appear when the nonlinear terms in the EL equations are included. 

Our presentation is organized as follows. In Sec. Ull we provide the physical background of 
the phenomena studied in the present work. In Sec. Illll we describe our model. Approximate, 
analytical solutions of the Euler- Lagrange equations for the order parameters, valid for 

> 1, are given and discussed in Sec. IIVI In Sec. |V] we obtain results for the effective 
potential between the confining surfaces. The quantitative comparison with the experimental 
data are described in Sec. |VT1 We discuss our results and conclude in Sec. IVIII 



II. BACKGROUND 



A wall of a container or a surface of a colloidal particle disturb the structure of the fluid 
in contact with them because of geometrical constraints on the positions of the molecules, 
and because the molecules interact with the matter of the wall rather than with the fluid 
molecules that are missing beyond the external surface. Structural changes are present for 
separations from the surface of the order of the bulk correlation length ^. In particular, near 
a surface preferentially adsorbing one component of the mixture the excess concentration of 
this component extends to distances ~ ^, and for T ^ Tc (hence ^ — 00) this phenomenon 
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is called critical adsorption 23 1. 



When a second, parallel wall at the separation L from the first one is present, the excess 



grand potential of the fluid confined in the slit has the form 24 1 



= uJexA = n + pAL = ((70 + 7l) + ^{L))A (1) 

where A is the area of each surface, p is the bulk pressure and 70, 7l are the surface tensions 
at the corresponding walls. The surface tension results from the particle-wall interactions, 
and from the modification of the structure of the fluid near the single surface in the semiin- 
finite system. The effective potential "^{L) reflects the mutual effect of both surfaces on the 
structure of the fluid. The structure of the fluid is influenced simultaneously by both walls 
if L ~ ^, therefore "^{L) vanishes for L ^ ^. Since the confined fluid tends to minimize 
the grand potential, "^{L) and — V\I^(L) act as an effective potential and an effective force 
between the confining surfaces respectively. Close to the critical point associated with either 
gas-liquid or demixing transition of the confined fluid ^ — )■ 00 and "^{L) acquires a universal 
contribution which becomes long-ranged at the critical point. This contribution to "^{L) is 
termed the critical Casimir potential and exhibits scaling described by a universal scaling 
function which is determined solely by the so-called universality class of the phase transition 
occurring in the bulk, the geometry, and the surface universality classes of the confining sur- 
faces. The range of "^{L) can be tuned by small temperature changes, because ^ = ^ot~'^, 
where t = {T — Tc)/Tc, is the critical temperature of the solvent, the critical exponent is 
u ~ 0.63 and the system-dependent parameter ^0 is of order of a few Angstroms. The theory 



of the thermodynamic Casimir force 
confinement 
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23 



based on the theory of critical phenomena in 
25j is well established. For the Ising universality class in a slit geometry, 
pertinent to the present study, "^{L) decays exponentially with the decay length equal to the 
bulk correlation length for distances L > C,- For the symmetrical (antisymmetrical) surfaces 
the potential is attractive (repulsive). In the case of a binary mixture, symmetrical (anti- 
symmetrical) surfaces have the same (opposite) adsorption preferences for the components 
of the binary mixture. Following a convention used commonly in the literature we denote 
by (+, +) (and, equivalently, (—,—)) the boundary conditions (BC) which reflect the fact 
that the two surfaces effectively attract the same component of a liquid mixture, whereas 
(+, — ) BC correspond to the case in which the two surfaces attract different components. 
According to the above discussion, one expects the effective critical Casimir interaction 
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to occur between a colloidal particle and a planar wall or between two colloidal particles 
immersed in a near-critical binary solvent. Often, in such systems also electrostatic interac- 
tions are present, e.g., in colloidal suspensions that are charge-stabilized. The charge at the 
colloidal particles or at the charged wall is screened by the counterions in the solvent. Ac- 
cordingly, the electrostatic interactions between two charged colloidal particles or between 
a colloidal particle and a charged wall become exponential functions of the distance and 
can compete with the critical Casimir forces. For instance, the effective interaction between 
charged planar surfaces decays as ±exp(— kL), where the repulsion (attraction) corresponds 
to the likely (oppositely) charged surfaces, and the dimensionless inverse Debye screening 
length is 



p* = pcO^ is the dimensionless number density of ions and a is the microscopic length 
unit (we shall choose for a the size of the solvent molecules). Moreover, this competition 
can become an interplay. In the present paper we consider charged surfaces immersed 
in a binary solvent. In such systems, if the solubility of ions in both components of the 
mixture is the same, then the distribution of charges is independent of the local solvent 
concentration, and also the concentration of the mixture is not affected by the presence of 
the ions. As a consequence, the presence of charges at the confining surfaces has no effect 
on the critical Casimir potential and, vice versa, the critical adsorption has no effect on 
the electrostatic interactions between the charged surfaces. Therefore ^{L) is just a sum 
of the critical Casimir and the electrostatic potentials. Usual salts, however, are soluble 
in water and insoluble in organic liquids. In such a case the critical adsorption of the 
component preferred by the wall and the distribution of charges may influence each other, 
and as a result may lead to a different form of ^{L). The room-temperature critical points 
are present in mixtures of water and organic liquids, therefore the question how the critical 
adsorption and the distribution of hydrophilic ions influence each other and modify ^{L) is 
of practical importance. 




(2) 
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III. GINZBURG-LANDAU THEORY 

In this section, following Ref. js] we briefly summarize the main steps in developing the 
Ginzburg- Landau theory both from a microscopic lattice gas model and from a continuum 
one. 



A. Derivation of the model 



In order to obtain the Ginzburg-Landau functional from a continuum microscopic model, 



one starts from the grand thermodynamic potential of the four- component mixture [24 1 

^ = UsR + Ue_i -TS- [ dr/i,p,(r), (3) 

Jv 

where Usr is the energy associated with the short-range (SR) vdW interactions, Uei is the 
electrostatic energy, S is the entropy, T is the temperature and /ij is the chemical potential 
of the i-th species. Local dimensionless number densities are denoted by p*(r) = pj(r)a^, 
where « = 1,2,3,4 for water, oil, -|- and — for ions, respectively. For a we have chosen the 
diameter of the organic molecules. In equilibrium, p*(r) correspond to the minimum of Q for 
given T, pj and the boundary conditions. For ionic species of the same valence, ps = /i4 = Pc 
because of the charge-neutrality condition. Integration (summation in the lattice version) in 
Eq. ([3]) is over the system volume V = AL, and summation convention for repeated indices 
is assumed in the whole paper. We assume the usual form of the internal energy Usr, 

Usr = AusR = j^^v dr'^p*(r)\/,,(r - v')g,,{r - r')p*(r') (4) 

rfrp*(r)r/(r). 



V 



where Vij and gij are the vdW interaction and the pair correlation function between the 
corresponding components respectively, and V^{y) is the sum of the direct wall-fluid po- 
tentials acting on the component i. In our model length is in a units, i.e. we consider 
dimensionless r* = r/a in (jll) and in the whole article. However, to simplify the notation we 
drop the asterisk for r as well as for the characteristic lengths (like see ([2])). It should 
be remembered that length is dimensionless. In the lattice model only nearest-neighbors 
interact, and the integration in Eq. (j4]) should be replaced by a summation. In contin- 
uum we assume that both the zeroth and the second moments, Vq^^ = J dTgij{r)Vij{r) and 
V^''*"'^ = ^ J drgij{r)Vij{r)r'^, respectively, are finite. 



Compressibility of the liquid can be neglected, so we assume 

Ei=iP* = 1- The three 
independent densities can be chosen as: a concentration of the solvent, 



Pi - P2> 



(5) 



a dimensionless density of the solute. 



Pc = Ps + P4 



(6) 



(subscript c from 'charge') and a dimensionless charge density. 



Ps - Pi- 



(7) 



Based on the experimental case where ions in the solution come from dissociation of a 
water, a similar chemical nature of the anion and the cation is assumed in Ref. 8|], and 
any difference between the interactions of the anion or the cation and any other species is 
neglected. In the case of salts insoluble in organic liquids the above assumption is not strictly 
valid, and should be considered as an approximation, whose validity should be verified at the 
later stage. This assumption distinguishes our analysis from Ref. 
consequence for the form of the short-range interaction energy. Namely, Usr expressed in 
terms of the new variables depends only on s and pc, and is independent of [8|], as can be 
verified by assuming Vi^s = Via in Eq. (jlj). 

The electrostatic energy in a slit with the surface charge cr(n) at the n-th wall (n = 0, L) 



101 ]. and has an important 



IS 



A 



Uel 



dz 



- — {Xjijf + ec, 



in 



(8) 



+ea{0)tlj{0) + ea{L)tlj{L), 



where e is the elementary charge, e is the dielectric constant of the solvent and the electro 
static potential ip satisfies the Poisson equation. 



+ echiz) 



0. 



(9) 



We neglect the dependence of e on the solvent concentration for two reasons. Firstly, we 
take into account that in the critical region the amplitude of the deviations from the average 
concentration is small an hence such a dependence leads to the higher order corrections to 
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the order parameter profiles, which we neglect (see Ref. 8| ). Secondly, as already men- 
tioned in the Introduction, the dielectrophoretic effects can be mimicked by an appropriate 
contributions to the surface fields. Accordingly, while comparing our results with the ex- 
perimental data we treat the surface fields as the fitting parameters. For an analysis of the 



dielectrophoretic effects see Refs 



• BE, 



11 



26|. 



The entropy S in the lattice model has the form of the ideal mixing entropy. Here we 
assume the same approximation. 



B. Separation of the charge-dependent and charge-independent parts of the grand 
potential 



The theory of critical phenomena was developed for uncharged systems, therefore we 
shall separate the part depending on the charge density from the remaining part of the 
grand potential. For the latter part we shall apply the Ginzburg-Landau description. 

In the new variables (Eqs. (jS]) and ([7])) S can be split into two terms, 

S = -ksA^^ / dzp*{z) In p*{z) = {sc[s,pc\ + Sei[pc, (p])A 
i=i -^0 



with 



sc[s,pc] = -ks I dz 
'o 



1 - Pc{z) + S(z) / 1 - Pc{z) + S{z) 



2 V 2 

- Pc{z) - S{z) f 1 - pc{z) - S{z) 



In 



-Pc{z) In 



Pc{z) 



and 



Sel[Pc 



-ks I dz 
'o 



Pc[z) + ^[Z) ( Pc[z) + (t)[z] 



Pc[Z) - (p[Zj_ ( Pc[Z) - (p[Z) 



Pc{z)\n 



Pc{z) 



(10) 



(11) 



(12) 



2 V 2 J ' ' \ 2 

We use the subscript 'eV for the quantities that are directly or indirectly associated with 
electrostatics and vanish for = 0, and the subscript 'C is from "Casimir". From the above 
properties it follows that the grand potential is a sum of the two terms 

Vt = {u]c[s,Pc\ +UJel[Pc,(p])A (13) 
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where 

Uel [pc ,4>\= Uel[4>\-TSel[PcA] (14) 

and 

i^c[s,Pc]=usr[s,Pc]-Tsc[s,Pc]- I dzpiPi{z) (15) 

Note that sc[s,p^ + A;Bln2 dzpc{z) (see Eq. flTTl) ) equals the entropy density of a three- 
component charge-neutral mixture with the solute density pc (i.e., Pz = P*i = Pc/'^) and the 
solvent concentration s in the case of close packing. Using this observation, one can interpret 
Eq. f|T5l) as the grand-potential density of such a three-component neutral mixture with pc 
replaced by pc + ksT ln2 (recall that pc = Ps = p^)- Thus, we have separated from the 
grand potential Q the contribution independent of the charge, which has this advantage that 
the methods developed for neutral near-critical systems can be directly applied to uc[s,pc\- 
Let us focus on Eq. (fT^ . which can be rewritten as 

UJel [pc, 0] = fel [pc, 0] " fcl [pc, 0] (16) 



where 



fel[pc, 0] = Uel[4>] + ksT / dz 

Jo 



p,{z))Hl-p,iz)) (17) 



Pcjz) + (j){z) / pcjz) +(j){z) \ ^ pcjz) - (f){z) fpc{z) - (j){z) 



This term alone describes the ions dissolved in a homogeneous solvent of the density 1 — 
Pc- However, the two contributions in Eq. f|T3|) . Uc and Uei, are coupled through pc- In 
equilibrium, pdz) corresponds to the minimum of Q (Eq. f ll3p ). therefore the critical Casimir 
and the electrostatic contributions to the effective potential are not independent, even in 
the case of identical chemical nature of the ions. This implies that the effective potential 
between the confining walls must differ from the sum of the critical Casimir potential and 
of the electrostatic potential in the case of a homogeneous solvent. 

C. Expansion of the functional 

For given T and pi the bulk equilibrium densities, s, pc and = 0, correspond to the 
minimum of the bulk part of fi, —pAL. Deviations of the fields s and pc from the bulk 
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equilibrium values are denoted by 

A{z) = s{z)-s (18) 
Mz)=Pciz)-p, (19) 

where z is the distance from the left wall. In equilibrium ^i{z), '&2{z) and (f){z) correspond 
to the minimum of Uex (see Eq. ([1])). Close to the critical point '(^i{z), 'd2{z) and (f){z) are 
small for z ^ C,, therefore the entropy can be Taylor expanded and the expansion can be 
truncated. From Eq. (fT2|) we have for fixed pc 



Sel[Pc + ^2,(P]= Sdh[(P] + As[?92, 



(20) 



where 







/ dz 
Jo 









(21) 



and 



As[i32,(p] = -kB [ dz 
Jo 



(22) 



. n>l m>l 

The above form follows from the fact that we have chosen to split the total entropy in such a 
way that Sei[pci 4>] vanishes for = 0. The coefficients a„ and an,m resulting from the Taylor 
expansion of Eq. f|T2|) are functions of pc- The above equations for fixed pc yield 



£e/[^2, 0] = UJel[Pc + ^^2, 0] = Cdh[(P] + AC[l92, 0] 



(23) 



where 



DH\ 



Uel\ 



ksT I dz 



(24) 



and 



AC[^2, 0] = -TAs[^2, 0] = -ksT [ dz 

Jo 



2p: 



(25) 



As already mentioned, Eqs. (I23ll- (l25l) with (jH]) and (19]) describe the ionic system with the 
charge density (p^z) and the total density of ions Pc + '&2{.z), placed between parallel charged 
walls. When the second term in Eq. fl23l) is neglected, no excess of the number density of 
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ions at the surfaces is obtained. It is the term ( 125|) of purely entropic origin that leads to 



the excess number density of ions near the surfaces |27l . |28| when the chemical nature (and 
hence the interactions with the wall) of the anion and the cation are the same. 

From the above considerations it follows that the excess grand potential can be split into 
three terms, 

Wex[^i,^2,0] ~ Cc[^uH + -CdhM + A£[i92,0], (26) 

where Cc[^i,'&2\ = ^c[pc + ^^2,5 + i^i] — wc[Pc, s]. Since Cc describes the near-critical two- 
component solvent with addition of one kind of neutral solute, it can be approximated by the 
Landau-type functional by using standard coarse-graining procedures. Close to the critical 
temperature 'di{z) and d2{z) vary on the length scale large compared to the molecular size, 
and 'di{z^),'d2{z') can be Taylor expanded about z' = z. From Eqs. (j4]) and ( ITT]) we thus 
obtain |8| 

Cc = C^c + ksT f dzY^Y. K.m^T{z¥2{z), (27) 

n m 

where in the summation in Eq. fl7r|) 2n + m > 3, hn^m are functions of pc and 

CcViM = \ ^'^<lA^,('A(flfi,[A + Vtf.(;)J„V!»,(z)| (28) 



where 



J°- = / drJijir) and Jij = | / dTJij{r)r^. —Jij{r) represents the vdW interactions for 'di and 
'dj, and can be obtained from the vdW contribution to Eq. ([3]) with the densities expressed 
in terms of the new variables (see ([5]) -([7])). We assume the same interaction ranges for all 
interacting pairs and postulate J^^ = 6Jij (recall that we consider dimensionless distance). 
Explicit expressions of C^j are given in Ref. and in Appendix A. Finally, 

hi{n) = hi{n) - Jij-djin), (30) 

where hi{n) are the surface fields describing direct interactions with the n-th wall. Jijdj{n) 
and the remaining surface terms in Eq. (|28|1 compensate for the interactions with the missing 
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fluid neighbors due to the presence of the wall; such interactions are present in the bulk term, 
but should be replaced by the interactions with the molecules of the wall js]. 

When the mixture phase separates, both the solvent concentration s and the density of 
the solute pc are different in the coexisting phases, because of a much bigger solubility of 
the solute in water. Likewise, for T close to both s and pc exhibit long-range critical 
fluctuations. Thus, in the Fourier representation the bulk part of Cq< can be written in the 
form 

A/:°,[i9i,^2] = AL0,[<fi,$2] = j rfk^[<li(-k)(7i(A;)l>i(k) + l>2(-k)C'2(A;)l>2(k)] (31) 

where Ci{k) and $i(k) are the eigenvalue and the eigenvector of Cij{k) = C°- + k'^Jij, 
respectively. The critical order parameter, $i(k), is associated with the eigenvalue Ci{k) 
that vanishes at Tc] C2{k) at Tc is positive and of the order of unity. The asymptotic decay 
of correlations in a real space is dominated by Ci(0) oc C,^^- In the critical region, the 
contribution from the noncritical OP $2(k) to the grand potential is much larger than the 
contribution from the critical OP. Accordingly, the probability of fluctuations associated 
with $2 is negligible compared to the probability of the fluctuations corresponding to the 
critical OP. This allows us to neglect the noncritical fluctuations and L^[$i,0] takes the 
usual form associated with the critical Casimir potential for the Ising universality class. 
From the computational point of view, in the present work it is more convenient to consider 
both fields, 'di and '1)2, instead of their linear combination $1. 



IV. APPROXIMATE SOLUTIONS OF THE EULER-LAGRANGE EQUATIONS 

In this section we derive the approximate EL equations for the functional ( l26l) and obtain 
approximate solutions for the solvent concentration, the solute density and the charge in a slit 
of width L ~ ^. In the one-phase region we neglect the second term in (1271) . and consider the 
lowest-order approximation which incorporates the coupling between the critical adsorption 
and the distribution of charges. 

The first, second and third terms on the RHS of Eq. (132!) are given by Eqs. (1281) . (12^ and 
fl25|) . respectively. We neglect the higher order terms in Eqs. flM|) and fl25|) . and in this 
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approximation the expansion of the functional is truncated at the third order term in the 
fields 'di and (j). 

The Euler-Lagrange equations, obtained by minimization of the approximate functional 
with respect to the fields di and </>, together with the Poisson equation Q, take a rather 



simple form 



dz"' 



M,^d,{z)+d4\z) (33) 



^ = "^^'^ + — d? — • 

In the above Mij = {J^^)ikC^j, where {J~^)ik is the (i,/c)-th element of the matrix inverse 
to the matrix Jij js], and (^1,^2) = —^^{{J~^)i2, {J~^)22^- The solutions must satisfy the 
charge neutrality condition, 

L 

dz(j){z) + (Jo + (Tl = 0, (35) 







and the boundary conditions for -^i: 

d^Az 



dz 
d^z) 

dz 

where 



2- 



.=0-^.(0) = i/.(0) (36) 



H^{n) = -{r%h,{n), (37) 

and h{n) is defined in Eq. (l30!) . For a hydrophilic (hydrophobic) wall Hi < {Hi > 0). 
Consistently with the approximate form of the functional fl32l) . the RHS of Eqs. fl33l) and 
(!34|) are truncated at the second order terms. 

A. Solutions of the linearized EL equations 

The linearized equations fl33l) and fl34|) for i)i and (j) are decoupled. The solutions take 
the well known forms 
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where we denote the ratio of the surface charges at the two surfaces by 

K = — (39) 

and use the superscript (1) to distinguish the solutions of the hnearized equations. The above 

charge profile obeys the charge neutrality condition fl5^ . The corresponding approximation 

for the electrostatic potential is (see hnearized Eq. and Eqs. 09]), (I2])) 

JcT 

= (40) 

Pc 

The excess concentration of the solvent and the excess number density of ions ( ITSl) in the 
critical region T — Tc (.^ — )■ oo) take the approximate form 

= toe-^/« + tie-(^-^)/«, (41) 

Because in the critical region the decay length associated with the larger eigenvalue 
(72(0) oc is negligible compared to ^, the terms oc e""^^ are subdominant and can be 
omitted for slits with L 3> a. From the boundary conditions we obtain the approximate 
expressions 

no - rioi ^ _ ^2g-2L/g (42) 



1 - CRte-^'^ 



to ~ toi 
ti ^ to 



1 - C2e-2^/« 
Rt - Ce-^/« 
1 - C2e-2i/« 



where 



and 



C^l^ (43) 



noi = =^c^oo -^2(0), toi ^a^oo -Hi{0), (44) 



(e+1^ 



_i72(L) i/i(L) 



if2(0)' //i(0)- 

We note that the decoupling of the fields '^2 and (p that occurs after linearization of EL 
equations is rather unphysical. Nonlinear terms are necessary in order to regain the right 
physics. 
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B. Leading-order corrections in the critical region 



In this section we determine the leading-order corrections to the solutions of the linearized 
EL equations fl55]) and fHT]) . In Ref. js] it was assumed that except from distances ~ a from 
each wall the dimensionless fields f = are all of the same order of magnitude, / = 0(z/), 
where z/ is a small parameter. and are proportional to Hi and a respectively, thus the 
analysis in Ref. is restricted to the surfaces with H^, a = 0{u). Under the above assump- 
tion analytical solution of the EL equations can be obtained by systematic approximations 
within a perturbation method. Since the RHS of Eqs. fl33|) and f p^ are truncated according 
to the truncation of the functional (see fl32l) ). in a consistent approximation the solutions 
should have the form -(^j = + df^ and 0j = (pf^ + The superscript (2) refers to 
the leading order correction terms (of order 0(z/2)), which satisfy the linear inhomogeneous 
equations 



^^.y-)(.) + l-^^^^^^^[;)5(-)). (46) 

In the above 0*^^^ and 'df'^ are given by Eqs. (138!) and (14T]) respectively. The boundary 
conditions are ddf\z) / dz\,=Q = 4'^\o), ddf\z) / dz\,=L = -^f\L) and j^; dz<P^^\z) = 0, 
because obey Eqs. (l36l) . and 0^^^ obeys the charge neutrality condition (l35l) . Note that 
because the Poisson equation is linear, from the above and ([2]) we obtain the leading-order 
correction to the electrostatic potential 



epc 



^2 9 _ ^(2) 
pc 



(47) 



Note that in a semiinfinite system the second term on the RHS in Eq. fl45|l decays as 
~ exp(— 2k2), and the second term on the RHS in Eq. (146!) decays as ~ exp{—Kz) exp(— 2;/^). 
Further approximations are possible when one of the two length scales, either the correlation 

nn 

or the screening length is much larger than the other length. Following Refs. [8|, |9[ 
we introduce the ratio between the correlation and the screening lengths, 

y = (48) 
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and focus on the two limiting cases: (i) y <^ 1, i.e., the Debye length is much larger 
than the correlation length, and (ii) y ^ 1, i.e., the Debye length is much smaller than 
the correlation length. The analysis of the limiting cases can be done with a reasonable 
effort. We should note that the experiments showing unusual attractive effective potential 
between the charged colloidal particle and the charged wall having the opposite adsorption 
preferences, were performed for ?/ < 1 [l^, whereas the experiments reported in Ref. 



concern the case y > 1. 



Is 



Let us first focus on the case (i), which was studied in Ref. [9|. For y <^ 1, from Eq. 
we have ~ exp(— ^ 1, and the second term on the RHS of Eq. ([31]) can be 

neglected. As a result we obtain that ~ 4'^^\ and "di satisfy Eqs. (1331) and (15^ . The 
solution of Eq. (155]) with (p ^ 0*^^^ yields a qualitative agreement with the experimental 
results for the effective potential obtained for a system with the Debye length larger than 



the correlation length 14 1. 



The case (ii) was studied in Ref. [8|] for a semiinfinite system. For y ^ I'm. the semiinfinite 
system we have (see Eq. (138]) ) {<t>''^\OY ~ exp(— 2y) ^ 1. Thus for z ^ ^ the second term 
on the RHS of Eqs. ( 133]) and ( 145]) is subdominant with respect to the second term on the 
RHS of Eq. din]) decaying as ~ exp(— kz) exp(— z/,^) exp(— y — 1). Therefore, in the 

asymptotic region of y ^ 1 we neglect the former but keep the latter. This means that we 
can approximate the concentration and the number density of ions by the solutions of the 
linearized EL equations di ~ but we cannot do it for the charge density 0. In physical 
terms the effect of the charge profile on the critical adsorption is negligible for ?/ ^ 1, because 
the neutralizing charge in the fluid is present at the distances from the surface z ~ ^ 
the charge distribution can be neglected on the same footing as the distribution of molecules 
at the distance from the wall. 

In this work we focus on the case of y > 1, and adopt the approximation valid in the 
asymptotic region y ^ 1. We neglect the effect of the charge distribution on -di and -^2, 
and obtain the leading-order correction to the charge profile from the approximate equation 
( 146]) . The solution of Eq. ( 146]) can be written in the form 

Ae""" + Ae-"(^-") + Aoe~"'(l - e~''^) + ALe^''"(l - e~'^''~'^'^] 



^ ^' Pc(l-e-^) 



(49) 
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z/L z/L 

FIG. 1: Charge density profile as a function of the scaled distance zjL. The solid line represents 
(Eq. ()38p . linearized DH theory result) and the dotted line is the approximate charge density, 
Eq. ()50p . with the effect of the critical adsorption included, = 5, kL = 5 , = \Rn\ = 1 (the 
same charge densities and the same or opposite adsorption preferences at both surfaces), and the 
excess number density of ions at the surface is noi = 0.5. Charge density is in (— kuo) units. Panels 
(a) and (b) correspond to the boundary conditions (— ,— ) and (— ,+) respectively. 

The coefficients are functions of no, ul and Ra-, and their explicit forms are given in 
Appendix B. The obtained approximation for the charge density 

(i){z) ^ <P^^\z) + 4)^^\z) (50) 

is shown in Figs. [I](a) and (b) for (— , — ) and (— , +) boundary conditions, respectively; 
results correspond to y = nL = 5. 

V. THE EFFECTIVE POTENTIAL 

The main goal of this work is a determination of the effective potential "^{L) between 
the surfaces that are both selective and charged. From Eq. fl32l) it immediately follows that 
the effective potential can be approximated by the sum of the Casimir and the electrostatic 
potentials only when the last term in fl32|) is neglected. When the chemical nature of the 
anion and the cation is the same, this term is of a purely entropic origin. 

When AC['d2,4'] in Eq. ([32]) is neglected, then the Casimir and the electrostatic terms 
are decoupled, and the concentration of the solvent and the solute density are obtained by 
the minimization of Cc['&i,^2], whereas the charge profile is obtained by the minimization 
of C£)h[(P] (with a simultaneous solution of the Poisson equation ([9])). 
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When the last term in Eq. (132|) is taken into account, it directly yields an extra contri- 
bution to the effective potential. What is important, this term depends on both, k, and C,, 
as well as on the surface charges and the surface fields. In addition, when AC['d2i4>] is in- 
cluded, then the electrostatic contribution is CdhI'P^^^ rather than C£)j{[(p^^'^] obtained 
in the absence of the critical adsorption. This term also depends on ^ and the surface fields 
through 0(2) (see (H9|l ). We stress again that for a homogeneous solvent A£['^2,0] leads to 
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28|. Neglecting this 



the excess number density of ions in the layer of thickness 1/{2k) 
term leads to an oversimplified theory already for a homogeneous solvent. 

In this section we determine the form of the potential "^{L) = Wex~7o~7L by substituting 
to Eq. (152]) the approximate forms of the fields -i^j ~ (Eq. (1411) ) and ~ 0(i) + 0^^) (see 
Eqs. fl55]) and fH^ ). The expression for ujex['&i,'^2, (t>] simplifies greatly when the fields 



§i,§2i<P satisfy the EL equations. For the Casimir part we obtain in our MP approximation 



(51) 



where Ac = —4:Hi{0)JijHj{L) (see Eq. fjTTl) ). The remaining contribution to Uex in the 
approximation consistent with Eq. (152]) can be written in the form (see 



£ei = *d(^) + 7eK0) + ML) ~ 41^ + 4?- 



(52) 



The leading order term (0(z/2)) is given by Eq. (IM]) with and approximated by the 
solutions 0(1) and -0(1) of the linearized equations. 



(1) 

el 



dz 



The well-known solutions are 



2pc ' 



,(1)2 



in 



V0(i)) +60(1)0.(1) 



+eao[^(i)(0) + i?.V^(i)(L)]. 



(53) 



and 



(1), 



n] 



DH 



Pc 



2pc 



1 f^L 
coth — 
V 2 



- 1 



(54) 



(55) 
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0.04- 
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□ The present theory 
— DH theory 




FIG. 2: The contribution to the effective potential per unit area associated with the presence of 
charges, ^ = ^^j^ + (Eqs. (j55p and (j58p ) for (— , — ) BC as a function of the scaled distance 
kL (dotted line). Solid line is the potential resulting from the linearized DH theory. The 



potential is in units of 



2pc 



kS, = 5, = Rn = I and noi = 0.5 (see Eqs.dMD and dS])). 



The leading-order correction term is of the order 0(z/^), and has the explicit form 



(2) 
el 



dz 



^0(1)2^(1) + ^0(1)0(2) _ + e(0(l)^(2) ^ ^(2)^(1)N 

2p2 47r 



(56) 



By using Eqs. f HOj) and f H7|) . integrating by parts and after some algebra we obtain 

kuT '•^ 



C 



(2) 



2p? 



+ e^(^)(L)|....o + ^'Ue^(^)(0)L ^^^^^(°) 



(57) 



where the first term equals A£['i92^V'^'']- The remaining terms come from + 
</)(^^]. (The full expression for Cf^ is given in Appendix C.) 

We neglect terms 0(exp(— 2L/^), exp(— 2fi:L)), and after subtracting the surface-tension con- 
tributions, we obtain the approximation 



2p: 



n2 



(58) 
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where the coefficients are (see Eqs. ( 139|) . (jH]) and (jUl 

K + Rnf y 



Ai(y) 



Ra \2y-l 



(59) 



A2{y) = {l + Rn) 



+ % 
2y+l 



(60) 



^3(y) 



-4(i + i?„) 



%2 



1 



(61) 



The above approximation is vahd for finite y; for ?/ — )■ oo it is not justified to neglect terms 
0(exp(— 2L/,^)), and the approximation (|32l) is oversimphfied. The case relevant for the 
experiments in Ref. 12, llSj, however, corresponds to 1 < y < 10. The effective potential 
+ (Eqs-(l55l) and (EHD) is shown in Fig. [2] for (-, -) BC. A similar electrostatic 
attraction between likely charged surfaces was found in Ref.[ll| in a nonlinear theory for 
large k^. 

The final approximate expression for the potential per unit surface area of the slit, in the 
region accessible in these experiments takes the form 



where (see Eq. (H4l) ) 



V^{y) = Acr' 



KaoaLUoi 



2pi 



My) 



(62) 



(63) 



T^2iy) 



Pc 



1 /I ^ ^ 



(64) 



T^3iy) 



2a 



Asiy). 



(65) 



In Eqs.dSSD-dSSD the excess number density of ions at one surface, noi = —H2{0), and the 
number density of ions in the bulk, pc, have the dimension of 1/ volume, surface number 
densities of elementary charges, <7o,(Tl have the dimension of 1/area (see Eq.([H])), and the 
inverse Debye and correlation lengths, k and 1/^ respectively, as well as the amplitude Ac, 
have the dimension of 1/ length. 
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The above MF result can be corrected, because as noted in Sec. IIIIBI the Casimir con- 
tribution to the potential can be considered separately, and the critical fluctuations can 
be incorporated in this part. The universal scaling function of the critical Casimir force 



has been obtained in Monte Carlo simulations 



21| . The amplitude Aq characterizing the 



long-distance decay of the potential has been extracted from the asymptotic behavior of 
this function for L/^ ^ 1 in Ref . [isl . In the case of the symmetrical BC ((+,+) or (-,- 
)) Ac = A^/C^ = — 1.51(2) /£, and in the case of the antisymmetrical BC ((+,-) or (-,+)) 



Ac = A./^ = 1.82(2)/^ [13|. 



Important consequence of the coupling between the critical adsorption and charge distri- 
bution is the dependence of the prefactors in Eq. fl62p on the ratio between the correlation 
and the screening lengths, y. 

VI. COMPARISON WITH THE EXPERIMENT 
A. Derjaguin approximation 

The theory developed in the previous sections concerns confining surfaces that are planar 



and parallel to each other, whereas the measurements in Refs. 



12 



13l ] with which we would 



like to compare our findings were performed for a planar substrate and a spherical colloidal 



particle. When the colloidal particle radius is much larger than the separation of its sur 



ace 



from the substrate, then the Derjaguin approximation can be applied, as in Refs.|12l. |13 |. 
The curved surface is approximated by a set of concentric circular rings of the infinitesimal 
area dS{6). The rings are parallel to the substrate and are at the normal distance L{6) = 
z + R{\ — cos^) (Fig. [3]). For each ring the excess grand potential per unit area is given 
in Eq. fl62l) . except that the surface charge ctr of the ring differs from op of the colloidal 
particle, and the relation between them is 

(Tr = ap/ cosO. (66) 

Consequently, the ratio between the surface charge at the substrate and at the ring is related 
to the corresponding ratio between the surface charge at the substrate and the particle by 

Ra[w/R] = Ru[w/P] ■ C0S6, (67) 
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w 



o 



R 



FIG. 3: Illustration of the Derjaguin approximation for the plate-sphere geometry, dp, an, 
are the charge density of the spherical particle, of the ring, and of the wall respectively, z is the 
minimal separation between the surface of the colloid and the planar wall. L{9) is the normal 
distance between the ring and the wall. 

where the symbols w,P,R denote the wall, the colloidal particle and the ring, respectively. 
The contribution of the ring to the potential between the substrate and the particle has the 
form 



where dS{6) is the area of the infinitesimal ring. Finally, the total potential "^{z) is obtained 



d^{z) = dS{9)'i'{L{9)) 



(68) 



by summing all the contributions d'^lz) of the circular rings up to the maximal angle 6m 

7r/2, 




(69) 



or 




(70) 
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B. Fitting 



In this section we shall compare the predictions of our theory with the experiments 



reported in Ref. 



12 



13| . In the experiment one surface was a charged surface of a particle, 



and the second surface was a flat, hkely charged substrate chemically treated to achieve a 
desired adsorption preference. 

Although the theory developed here is of the mean-field type, the Renormalization Group 
(RG) results can be applied to the Casimir part according to the discussion in Sec. IIII B[ We 
shall assume that the general form of the potential, Eq. fl62|) . is a fair approximation, except 
that the correlation length should have the correct temperature dependence, ^ = ^o''""'^; i-e., 
with V taking the three-dimensional value 0.63 of the Ising universality class. Moreover, we 
shall assume that the Casimir amplitude Aq is given by the proper universal form associated 
with the Ising universality class. 

We shall compare our predictions with the experiment for all four pairs of the boundary 
conditions: (+,+), (—,—), (+, — ) and (—,+), where (+) denotes a hydrophobic and (— ) 
denotes a hydrophilic surface; the left and the right symbol in the pair refer to the particle 
and the substrate respectively. Unfortunately, neither the charge density = oi^ nor the 
surface fields Hi{L),i = 1,2 could be measured experimentally. Two kinds of colloidal 
particles were used: a hydrophilic with the unknown charge density and the radius R = 
1200(nm), and a hydrophobic with the radius R = 1850(?im). The surface fields Hi{0),i = 
1,2 characterizing the colloidal particles are also unknown. According to the experimental 
conditions we assume that the surface charge ap = ctq and the surface fields for the colloidal 
particles of the same type are fixed. We thus impose strict constraints on H2{0) and ctq 
to take on the same values for (+, +) and (+, — ) BC, and the same values (but of course 
different than in the previous case) for (— , — ) and (— ,+) BC. The same constrains are 
imposed for the parameters that describe the fiat surfaces, i.e., we require that cr^ = croi?o- 
and H2{L) = H2{0)Rn are the same for (— , — ) and (+, — ) BC, and likewise the same for 
(+,+) and (-,+) BC. 



In experiments of Refs. 12|, ll3| , ions in the solution were present due to water dissociation 



in a salt free water-lutidine mixture. For this mixture, according to Ref. 29|] the density of 
(monovalent) ions is about pc — 1.08 • 10~'^mo///. We use this value, although it appears to 
be a rather rough estimate, and consider fitting parameter. 
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FIG. 4: (Color online) The effective potential between a wall and a colloidal particle of the radius 
R = 1 200 (nm) immersed in a water-lutidine mixture as a function of the distance z for various 

n 

temperatures T [13[. The data refer to (— , — ) BC corresponding to the case where both the colloidal 
particle and the substrate are hydrophilic. Here Tc is the critical temperature of the mixture. The 
solid lines are the theoretical predictions (Eqs. (j70p and (j62p ) as explained in the main text. The 
parameters obtained from the fitting are shown in Table HI 

We take into account that the wall-particle distance z was determined in experiments up 
to zq = ±30(nm); we assume that for the specific boundary condition (the same series of 
measurements) the shift between the actual and measured distance is fixed; the shift can 
differ from one series of measurements to another. In the fittings, we have tried to keep 
the same value for for all sets of boundary conditions. The amplitudes Ac = 2tiA±R 



have been taken from Ref. [I3j; Aj^ and A_ are the amplitudes governing the asymptotic 
decay {L/^ ^ 1, r > 0) of the universal scaling functions of the critical Casimir force for 
symmetrical and antisymmetrical BC in a slit, respectively; for the Ising universality class 
in (i = 3, the Casimir scaling functions in a slit geometry were obtained by MC simulation 
method in Ref. 211]. 



In Fig. m we show the comparison between our theoretical predictions, Eq. fITOl) with 
\1/ given by Eq. (!62l) (solid lines), and the experimental data of Ref. [3| for (— , — ) BC. 
The obtained fit parameters are given in Table [B According to the table, the values of the 
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Tc-T 


(,{nm) 


Ac{nm) 


cro(nm) ^ 


aL{nm) ^ 


if2(0)(nm)-3 


H2{L){nm)-^ 


K.{nm) ^ 


zq (nm) 


0.14 


25.9 


-11379 


1.28 


0.064 


-0.01 


-0.0002 


0.0872 


-30 


0.16 


24.3 


-11379 


1.28 


0.064 


-0.01 


-0.0002 


0.09 


-30 


0.18 


23.3 


-11379 


1.28 


0.064 


-0.01 


-0.0002 


0.091 


-30 


0.19 


22.85 


-11379 


1.28 


0.064 


-0.01 


-0.0002 


0.0918 


-30 


0.21 


21.24 


-11379 


1.28 


0.064 


-0.01 


-0.0002 


0.094 


-30 


0.30 


16.3 


-11379 


1.28 


0.064 


-0.01 


-0.0002 


0.097 


-30 



TABLE I: Fit parameters for the effective potential given by Eqs. (170p and (162p for (— ,— ) BC 
where the colloidal particle with the radius R = 1200(nm) and the substrate are both hydrophilic. 
The amplitude Ac for this system is taken from Ref. isj. ctq and ai denote the surface charge at 
the particle and at the substrate respectively in units of elementary charge e. -^2(0) and H2{L) 
denote respectively the dimensionless effective potential per unit volume between the particle and 
ions, and the flat substrate and ions (see Eqs. ()37p and ()44p ). zq is the experimental error in the 
measured distance between the substrate and the particle. and k are the correlation and the 
inverse Debye-length respectively. See the main text for more details. 



correlation length 



fit 



rpfit 



-0.63 



(71) 



with ,^Q** = 0.21 ± 0.004(nm) and the Debye screening length = 10.9 ± 0.6(nm) are 
both in the range of the experimental results. The best fit is obtained for a zero shift 
in the critical temperature, T/** = but taking into account up to 5mK inaccuracy in 
T itself. The charge density of the colloid (in units of the elementary charge e) obtained 
from the fit is ctq = 1.28(r2m)~^, which agrees nicely with the value given in experiments 



of Refs. 



29 



30|, and is compatible with the observation that the highly charged colloids 
(> 0.24(nm)~^) preferentially adsorb water while the colloids with a smaller amount of 
charge prefer lutidine. 

Figure \5\ shows the experimental data (symbols) and the theoretical curves (solid lines) 
for the case of (-I-, — ) BC where the colloidal particle is hydrophobic whereas the wall is 
hydrophilic. The obtained fit parameters are given in Table Ull According to the table, 
the amplitude of the correlation length and the Debye screening length are estimated as 
■Co** = 0-21 ± O.OOl(nm) and = lO(nm), respectively. H2{L) and ai are the same as for 
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(+,-) BC 




0.4 

z(|am) 



FIG. 5: (Color online) The same as in Fig. Hlbut for (+, — ) BC corresponding to the hydrophobic 
colloidal particle and the hydrophilic wall. The radius of the colloidal particle is R = 1850(nm). 
The obtained parameters from the fitting are shown in Table HIl 

the (— , — ) BC. In this case, the best fit is obtained with a shift in the critical temperature, 
AT/^* = |Tc - T/**| = 223mK and allowing up to 5mK inaccuracy in T itself. The charge 



density of the colloid (in units of e) obtained from fitting is ctq = 0.7 ■ 10 



[nm] 



which is 



in agreement with Refs. 29[ and 31[ reporting the values for the surface charge densities of 
silica and polystyrene spheres in water. 



Similarly, for (— , +) BC the comparison of the experimental data of ISj and our theoret- 
ical predictions for the effective potential are shown in Fig. [61 The obtained fit parameters 
are given in Table UTTl For this case, the estimates for the values of the amplitude of the corre- 
lation length and the Debye screening length are ^q** = 0.21 ±0.003(nm) and = lO(nm), 
respectively. As for the (— , — ) BC, the best fit is obtained for considering no shift in the 
critical temperature, T/** = T^, and up to 5mK inaccuracy in T itself. The charge densi- 
ties and the surface fields of the colloidal particle and the substrate are consistent with the 
results of fitting for (— , — ) BC. 

For the experimental data corresponding to the (-I-, -|-) BC, i.e., where both the colloidal 
particle and the wall are hydrophobic, we find that 1.1 < < 2.1 (see Table HVl) . which 
means that the approximation >> 1 under which we have obtained Eq. ( l62l) is not 
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Tc-T 


(,{nm) 


Ac{nm) 


cjo(nm) ^ 


aL{nm) ^ 


if2(0)(nm)-3 


H2{L){nm)-^ 


K.{nm) ^ 


zq (nm) 


0.25 


85.77 


21144 


0.0007 


0.064 


0.001 


-0.0002 


0.1 


-18 


0.28 


49.5 


21144 


0.0007 


0.064 


0.001 


-0.0002 


0.1 


-18 


0.30 


39.8 


21144 


0.0007 


0.064 


0.001 


-0.0002 


0.1 


-18 


0.32 


34.8 


21144 


0.0007 


0.064 


0.001 


-0.0002 


0.1 


-18 


0.34 


29.4 


21144 


0.0007 


0.064 


0.001 


-0.0002 


0.1 


-18 


0.43 


20.5 


21144 


0.0007 


0.064 


0.001 


-0.0002 


0.1 


-18 



TABLE II: Fit parameters for the effective potential given in Eqs. (170p and (j62p for — ) BC 
corresponding to the hydrophobic cohoidal particle (with the radius R = 1850(n?7i)) and the 
hydrophilic wall. The amplitude Ac for this system is taken from Ref. For the remaining 

parameters see the caption of Table [H 
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FIG. 6: (Color online) The same as in Fig. [3] but for (— ,-|-) BC corresponding to the hydrophilic 
colloidal particle and the hydrophobic wall. The radius of the colloidal particle is i? = 1200(?im). 

strictly valid. In this case, the terms which we have ignored in Eqs. f l45p and f l46p (the terms 
which decay as exp(— 2k2;)) should be kept for a better comparison with the experiment. 
Nevertheless, we have performed fitting employing the relatively simple approximate form 
( 162|) of an effective potential. Because the neglected terms play a more significant role for 
small distances, in Fig. [7] we report the comparison with the experimental data only for 
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Tc-T 




Ac{nm) 


(To (nm) ^ 


aL{nm) ^ 


if2(0)(nm)-3 


H2{L){nm)-^ 


Hi{nm) ^ 


zq {nm) 


0.04 


55.07 


13715 


1.28 


0.1344 


-0.01 


0.0001 


0.1 


30 


0.05 


50 


13715 


1.28 


0.1344 


-0.01 


0.0001 


0.1 


30 


0.07 


42.5 


13715 


1.28 


0.1344 


-0.01 


0.0001 


0.1 


30 


0.09 


34.64 


13715 


1.28 


0.1344 


-0.01 


0.0001 


0.1 


30 


0.12 


29.64 


13715 


1.28 


0.1344 


-0.01 


0.0001 


0.1 


30 


0.18 


23.50 


13715 


1.28 


0.1344 


-0.01 


0.0001 


0.1 


30 



TABLE III: Fit parameters for the effective potential given by Eqs. (j70p and (j62p for (— ,-|-) 
BC, where the colloidal particle with radius R = 1200(rim) is hydrophilic while the substrate is 



131 ]. For more information see 



hydrophobic. The amplitude Ac for this system is taken from Ref. 
the caption of table I. 

distances larger than z = 80 (nm). The obtained fit parameters are given in Table IIV[ 
According to the table our estimate for the amplitude of the correlation length is ^q** ~ 
0.21 ± 0.018(nm) whereas for the Debye screening length we have k,~^ = 14.35 ± 0.6(?7,m). 
The smaller value obtained for n for this BC indicates, according to Eq. ([2]), that the ion 
density pc is smaller compared to the other BC. Since the experiment was performed for 

out of the range of validity of our approximate result, further studies for ~ 1 are 
required to verify whether the neglected terms in the potential would lead to the fit with 
~ lOnm. For the present case, the best fit is obtained with a shift in the critical 
temperature AT/** = |Tc — T/**| = 63mK and up to IbmK inaccuracy in T itself. Again, 
the charge densities and surface fields of the colloidal particle and the substrate are consistent 
with the fitting parameters obtained for (— , -|-) and — ) boundary conditions respectively. 

We find that the best fit for for all BC is £n = 0.21(?7,m). This value is in a very good 



agreement with several experimental results 
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. The correlation length ^ obtained from 



fits given in Tabs Hl lIVI are compared with the expected behavior 
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FIG. 7: (Color online) The same as in Fig. S] but for (+, +) BC corresponding to hydrophobic 
both the colloidal particle and the wall. The radius of the colloidal particle is R = 1850(nm). The 
obtained parameters from the fitting are shown in Table IIVI 



Tc-T 




Ac{nm) 


cJo(n?n) ^ 


aL{nm) ^ 


i72(0)(nm)-3 


i72(i)(nm)~3 


K{nm) ^ 


Zo{nin) 


0.08 


30.1 


-17543 


0.0007 


0.1344 


0.001 


0.0001 


0.07 


-15 


0.1 


25.93 


-17543 


0.0007 


0.1344 


0.001 


0.0001 


0.0675 


-15 


0.11 


23.8 


-17543 


0.0007 


0.1344 


0.001 


0.0001 


0.0669 


-15 


0.12 


21 


-17543 


0.0007 


0.1344 


0.001 


0.0001 


0.069 


-15 


0.13 


19.3 


-17543 


0.0007 


0.1344 


0.001 


0.0001 


0.0708 


-15 


0.2 


15.8 


-17543 


0.0007 


0.1344 


0.001 


0.0001 


0.0727 


-15 



TABLE IV: Fit parameters for the effective potential given given by Eqs. (I70p and (j62p for (-I-, -|-) 
BC, where the colloidal particle with radius R = 1850(nm) and the substrate are both hydrophobic. 
Only the data larger than z = 80{nm) has been considered. The amplitude Ac for this system is 
taken from Ref. |13i]. For the remaining parameters see the caption of Table I. 

VII. SUMMARY AND DISCUSSION 

In this paper we have studied a mutual effect of the critical adsorption and the distribution 
of ions on the effective potential between charged and selective surfaces confining a near- 
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FIG. 8: The solid line shows the correlation length ^, Eq. (j72p with the amplitude = 0.21 (njn) as 
a function of the deviation AT = Tc — T from the critical temperature Tc- Symbols correspond to 
the best fits of the effective potential given by Eqs. ([70]) and ([62]) (Table Hl UVl) . Recall that the case 
of (+,+) BC is out of range of validity of our approximation, and good agreement is not expected 
for this case. 

critical binary mixture with ions. We have employed a Ginzburg-Landau-like theory which 
can be derived either from a lattice gas model jsl or from a simple density functional 
theory jol for a four component mixture. We assumed the same chemical nature of the 
anion and the cation, and a much bigger solubility of the ions in a one component of a 
binary solvent than in the other. Such conditions are met, e.g., for aqueous solutions in 
which the ions come from a dissociation of a water. We have focused on the vicinity of 
the critical point of the solvent, where the correlation length for fluctuations of the solvent 
concentration, ^, is much bigger than the screening length I/k. 

We have shown that when the chemical nature of the ions is the same, the excess grand 
potential of a system confined between two parallel walls can be split into two parts, a part 
which is independent of the charge distribution, and a part which is independent of the 
solvent concentration. The first contribution describes a near-critical binary solvent with 
a neutral solute (uncharged ions ) comprised of single species with a preferential solubility 
in water. This part yields the critical Casimir potential. The second contribution to the 
excess grand potential is associated with the charge distribution, and has the form known 
from the DH theory. Both contributions, the critical Casimir and the DH, depend on the 
number density of ions. The equilibrium form of the number density of ions, corresponding 
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to the minimum of the grand potential, is different from that which arises from the critical 
Casimir part alone, and from the DH part alone. Thus, the effective potential between the 
confining walls differs from a sum of the Casimir potential in the uncharged system and the 
DH theory prediction for the electrostatic potential for ions in a homogeneous solvent. 

We have shown that for ^ 1 the effect of the critical adsorption on the charge distri- 
bution dominates, and the effect of charges on the solvent concentration can be neglected 
8| . This is because the screening length is much shorter than the correlation length, and the 
charges that are present at distances from the wall much smaller than ^ can be neglected, 
like the other molecular details. Because of the preferential solubility in water, the excess 
number density of ions decays in the same fashion as the excess solvent concentration, i.e. 
~ exp(— 2;/,^). The number density of ions determined by the critical adsorption is an input 
in calculations of the electrostatic contribution to the effective potential. It influences the 
charge distribution, and in addition changes the entropic contribution in the DH part of the 
grand potential. As a result we obtain terms which were absent in the standard DH potential 
in the case of the homogeneous solvent. The dominant additional terms are ~ exp[—L/^) 
and ~ exp(— kL). 

Note that since the critical Casimir contribution to the excess grand potential has the 
same form as in the uncharged critical system and since for ^ 1 it can be considered 
independently of the remaining contribution to the excess grand potential, the RG theory 
can be applied to this part. Thanks to the above separation and thanks to the analytical 
solution, we have been able to incorporate in our theory the RG results for the critical 
Casimir part. The remaining electrostatic part was obtained on the MF level, with correct 
form of the number density of ions resulting from the critical adsorption. The key result of 
the RG theory is the universal form of the critical Casimir potential, depending only on the 
boundary conditions as discussed in the Introduction. The molecular details, described in 
this theory by the vdW interaction potentials, influence only the nonuniversal properties, 
in particular the amplitude .^0 of the correlation function and the amplitude of the excess 
solvent concentration and the density of ions. 

The dominant term in the electrostatic contribution to the excess grand potential decays 
in the same way as the critical Casimir potential. This electrostatic contribution is nonuni- 
versal. Even though the critical Casimir part of the excess grand potential should exhibit 
the universal behavior, it can be obtained from experiments only when the electrostatic 
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contribution is subtracted. 

In order to verify the theory, we have fitted our predictions to the experimental results 



12 



13| . Two different substrates and two different particles were used in experiments to yield 



4 combinations of the boundary conditions. We have kept the same values of parameters 
characterizing the same surface. This requirement has provided us a constraint on the fitting 
parameters for the surface charge and for the surface fields. We have used the Derjaguin 
approximation to take into account the curvature of the surface of the colloidal particle. 
We have obtained a good quantitative agreement for a large range of distances (Figs. H]- 
in three cases, and a less good agreement (Fig. [7]) in the fourth case, which is at the 
limits of applicability of the approximations we have made. In our fitting, the amplitude 
^0 = 0.21nm has the same value for all the considered cases, and this value is in a very good 
agreement with various experimental estimates 32|-l36l|. The best fit is obtained for the 



correlation length that agrees very well with the expected behavior (see Fig. |8]) (except from 
the boundary conditions (+,+), where the agreement is less good). Also the fitted value of 
the critical temperature was precisely equal to the experimental value for (— , — ) and (— , +) 
BC, with small shifts, AT^ = 63mK and ATc = 223mK, for the (+, +) and (+, -) BC, 
respectively. 

As already rnentioned in the Introduction, the attempt to fit experimental results reported 



in Ref. 



12 



13| to the sum of the critical Casimir and the electrostatic potentials failed. 



Moreover, in order to fit the experimental results to the universal Casimir potential for 
large separations, a different value of had to be chosen for each BC, from 0.1 7nm to 
0.26nm for (— , — ) and (— , +) BC respectively, and the (+, +) BC could not be fitted to the 
Casimir potential alone. A smaller value of leads to a larger prefactor Aq/C, multiplying 
the exponential decay exp{—L/^) of the potential. In the present theory the prefactor of 
the decay exp{—L/^) contains the electrostatic contribution in addition to the universal 
Casimir amplitude. Therefore both, the present theory with = 0.21nm and the pure 
critical Casimir potential with = O.lTram can yield a good fit at large distances for (— , — ) 
BC. However, unlike in the present approach, the other BC could not be fitted by the critical 
Casimir potential alone with the same value of .^o; a-nd the relative difference between the 
fitted amplitudes was as large as (0.26 — 0.17)/0.17 > 50%. We cannot find explanation 
for such large differences in the amplitudes for essentially the same mixtures. Boundary 
conditions should not have any effect on the bulk properties. Moreover, in Ref. 12|, |l3| only 
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distances much larger from the position of the potential minimum could be fitted, whereas 
the present theory yields a good quantitative agreement for a wide range of distances. The 
agreement is obtained for all the measured systems for the parameters that agree very well 
with the experimental data, and if precise experimental data are absent, are of a correct 
order of magnitude. There are no data for the surface fields of the four surfaces. The 
four free parameters, however, satisfy all the constraints of consistency for the four pairs of 
surfaces (24 curves in 4 series of measurements). 

We conclude by stressing that a very important advantage of the analytical expression 
(Eqs.f lH^ -f lB^ ) is the possibility of designing the effective potential of a desired form by 
adjusting the surface charges and/or hydrophilicity (or hydrophobicity) of the surfaces, or 
the amount of ions in the solution. Eq.f lB^ may be a very useful tool in guiding future 
experimental studies. 

Acknowledgments 

We have greatly benefited from discussions with Siegfried Dietrich, Ursula Nellen, Marcus 
Bier and especially Laurent Helden, who explained us the details of experiments and com- 
mented on the results of fitting. Faezeh Pousaneh would like to thank Prof. Dietrich and his 
group for hospitality during her stay in Stuttgart, where a part of this work was done. The 
work of Faezeh Pousaneh was realized within the International PhD Projects Programme of 
the Foundation for Polish Science, cofinanced from European Regional Development Fund 
within Innovative Economy Operational Programme " Grants for innovation" . 

VIII. APPENDIX A. EXPLICIT EXPRESSIONS FOR C% 

The coefficients in Eq. (I29l) take the exphcit forms 
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IX. APPENDIX B. EXPLICIT EXPRESSIONS FOR THE COEFFICIENTS IN 
EQ. (1491) 



The coefficients in Eq. f l49p take the exphcit forms 

\y+lj [ l-e-'- J \y-lj 



and 



where 



At 



A 



LL 



y+l 



l-y- 



-kL 



l-e- 



-kL 



A 



LO 



Aoo = no 



{y + l? 

(21/ + 1)' 



AoL = -riL 



(2y - 1) 



Alo = -no 



Ra{y-i] 

(22/ -1) 



All = nL 



R.{y + 1)' 



l-y- 



1 - e-^/« 



1 - e-*" 



kL 



l-y- 



1 - e-^/« 



l-e- 



-kL 



(76) 



(77) 



(78) 



(79) 



(80) 



(81) 



(2y + l) 

Note that each coefficient (178l) - fl8Tl) diverges for y = ^ oo. However, the term 
in (1491) remains finite, because when k — oo then ye '^^ — 0, and when 
,^ — oo then y(l — e~^^^) ~ yz/^ = kz. Likewise, the whole correction to the charge proffie, 
Eq. (09]), is finite. 



X. APPENDIX C. EXPLICIT FORM OF THE LEADING-ORDER CORREC- 
TION TERM EQ. dH]) 

'(2) • 



Full expression for the in Eq. (157]) has following form 
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and the expressions for and nj;, are given in Eq. 
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